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Abstract 

We study the states localized near an interface between conducting and insulating bilayer 
graphene (BLG) and show that they have highly unusual properties that have no analog in con- 
ventional systems. Moreover, the states belonging to the two independent valleys in the Brillouin 
zone of BLG show contrasting properties that allow an easy experimental realization of various 
valley based functionalities desired in valleytronics without requiring any sophisticated techniques. 
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Bilayer graphene — two layers of carbon atoms with the same intralayer arrangement 
of atoms and relative stacking of the layers as in natural graphite, is focus of intensive 
experimental and theoretical research these days due to its unusual physical propertiesjl-8]. 
There are two independent valleys in Brillouin zone of BLG and the low energy elementary 
excitations in both valleys posses a pseudospin and have pseudospinor wavef unctions. The 
valley degree of freedom can be used for controlling an electronic device similar to the 
spin or charge and there are many proposals for achieving basic functionalities desired in 
this field, named valleytronics, using graphene systems {9 -14]. In refjioj], A. Rycerz et. al. 
proposed a valley filter that is a strip of monolayer graphene containing a narrow constriction 
with zigzag edges. However, experimental realization of such a device is challenging, if not 
impossible. Taking advantage of the effects of trigonal warping at high carrier densities in 



monolayer graphene, J. L. Garcia- Pomar et. al.[ll| proposed a device which again requires 
zigzag edges in addition to manipulation of beams of particles— again a great challenge to 
the experiment. In bilayer graphene, A. S. Moskalenko and J. Berakdar 12J showed that 
intense shaped light pulse can induce valley polarized currents and, D. S. L. Abergel and 
Tapash Chakraborty 13( proposed a device for generating valley polarized currents where 
in the presence of an intense terahertz light source in a region of bilayer graphene with a 
finite band gap the dynamically induced states can be made to exist only in a given valley 
by tuning different parameters. Although, these are relatively easy for experimental tests, 
they are hard to implement in practical devices. H. Schomerus' proposal 14j uses the fact 
that the angular dependence of the transmission probabilities of the particles of the two 
valleys from one region to another differ when the two regions have opposite pseudospin 
polarization. It demands a relatively sharp fermi surface and again controlling beams of 
particles and their propagation angles. Here we show that states that are localized near an 
interface of zero band gap and finite band gap BLG have highly unusual properties that 
provide us with a conceptually extremely simple way of obtaining valley polarization and 
related basic functionalities, and eliminate almost all experimental challenges that stand in 
the way of valleytronic devices. 

In a clean zero band gap or conducting BLG at any energy there are two pseudospinor 
plane wave states and two pseudospinor evanescent wave states. The latter have appreciable 
magnitudes only near interfaces between regions of different potentials and play an important 
role in transport properties. Similar is the case of a clean finite band gap or insulating BLG 
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where the evanescent waves are localized near interfaces between regions of different band 
gaps or potentials. For the first case they are only present at energies which are outside the 
band gap of at least one side. For energies outside band gap these states have the usual 
properties, however, for energies inside the band gap, their wavevectors become complex so 
the spatial parts of their wavefunctions contain plane waves as well as evanescent ones. This 
is unusual for any wave in a pristine and dissipation-free system and here it is more unusual 
as it apparently shows that number of particles is not conserved. In search of resolution to 
this issue we derive the continuity equation associated with the Hamiltonian of the system 
and check the current these states carry along the decay direction. It comes out to be zero 
just like the usual evanescent waves which is quite interesting. However, these states show a 
more interesting behavior in the transverse direction: the states belonging to the two valleys 
carry currents in opposite and fixed directions, independent of the sign of the wavevector 
component along this direction. The latter property is in stark contrast to common systems 
described by the usual Schrodinger wave equation. For an interface of a conducting and 
insulating BLG, due to the continuity of probability current density, these states modifies 
the behavior of the states localized on the conducting side as well. Thus the particles of a 
given valley in conducting region at any energy inside the band gap of the insulating BLG 
coming towards the interface at any angle turn towards the same lateral side close to the 
interface. Thus without worrying about the propagation angles of the quasiparticles or the 
finite smearing of the fermi surface due to temperature or any other intravalley scatterings, 
we can easily separate particles of the two valleys and also perform other valley based 
operations just by using a simple junction of conducting and insulating BLG either in above 
mentioned geometry or by applying a voltage difference along it as proposed in ref Q for the 
topologically confined zero energy chiral modes. 

The low energy electronics properties of the BLG system lying in the xy-plane can be 
approximately described by effective Hamiltonians H ± for the two valleys K(+) and K'(-) 
3, ISj which along with their pseudospinor eigenf unctions \l/ ± (a;,?/) can be written as 

H± =\ V U ± ^W(^)H r !l a) 

Where the components an d of 1 & ± (x,y) are envelope functions that give the 
amplitudes on the two different layers on sites A and B that do not lie directly above or 
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below sites on the other layer. V± = ^-{d x ± idy) 2 where d x>y are differential operators, 
the two layers have on-site energies equal to ±U induced by applied electrostatic gates, and 
m is the effective mass of quasiparticles when U = 0. When U ^ inversion symmetry 
between the two layers is broken and a band gap equal to 2\U\ is produced. Dispersion 



for both valleys is the same and reads: J5(k) = ±yU 2 + (^-) 2 where ± refer to electron 
and hole bands here. Suppose an interface between conducting region (CR, x < 0) and 
insulating region (IR, x > 0) of BLG at x = then U — for x < and we have continuum 
of propagating and evanescent states in CR while in IR for energies inside the gap only 
the states localized near the interface are present. The states in CR match those in the 
IR with the same energy E and, due to the invariance along the interface, with the same 
y-component of wavevector k y to ensure the continuity of the probability density and the 
probability flux density. Thus we can use E and k y to label all the states. In IR, for a given E 
and k y , the x-components of wavevectors ±k x are complex quantities given by k x = ±k r ±iki 
where k r = rjsma and ki = r/cosa, and we defined r\ = (5 2 + ky) 1 / 4 , a = | tan~ 1 (5//c 2 ) and 
5 = \/U 2 — E 2 . For simplicity we absorbed ^ in energy terms and we will do so hereafter. 
This results in decaying plane waves. The wavefunctions of the states belonging to the two 
valleys are ^/f(x,y) = ^/f(x)e tkyV where, ignoring the solutions that diverge with x, ^ff(x) 
are given by 




E-U \ 

q±(x) = af\ p n± _ ] e ikrX - kiX + 4 ( fl± n± 2 . 1 e - ikrX ~ kiX (2) 

Here af and df are complex coefficients determined by matching conditions, R± = 
k r 2 + ± ky) 2 , Q± = y/(E - U) 2 + H 2 ± and 0± = tan-^^^). The spatial parts of the 
above wavefunctions apparently show currents along ±x with magnitudes decreasing with x, 
however, as we will shortly confirm, the relative phases ±26 ± between upper and lower com- 
ponents of the pseudospinors contain the information needed to produce the proper results. 
We derive expressions for probability current densities J 1 * 1 for the two valleys associated with 
the Hamiltonians H in equation!! they read: J = (J^,Jy) where and are given 
by 

J± = ^{Im(-^a x d y V ± &a y d x *)} 
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Here a X:V are Pauli's matrices. It is also worth mentioning that these expressions show that 
unlike common systems the continuity of flux density is ensured with that of wavefunction 
and its first derivative in spite of the fact that in case of finite band gap BLG we have 
energy dependant effective mass as is clear from the dispersion relation. Further, the mixed 
x, y terms shows some intimate link between the two orthogonal directions. Using these 
expressions and defining (3 = tan -1 ^) the angle that q = (k r , k y ) makes with the x-axis, it 
is straightforward to show that the currents Jf = (4; 4y) carried by \4 (x) are 

J± = -6^(14 1 2 - |4| 2 )gcos(2fl ± =F (3)e' 2k ^ 

4 = C^KI 2 + \df\ 2 )qsm(26 ± T 0) + 2k y Re(44*e 2 ^ ±+fc -)) }e~ 2 ^ 

where C 1 * 1 = A ^ U ~^} R± , q = |q| and * shows complex conjugation. Note that the cross 
terms in J^ x vanish. Using a few common trigonometric identities we get 26 ± =p P = vr/2 

independent of any parameter so both terms in J^ a vanish to give 

4 a = o 

Thus the relative phases between components of pseudospinors count for the presence of 
plane waves and we obtain sensible results. A more interesting result is the currents along 
y-direction: 



J% = ^{Td^l 2 + |4| 2 )g + 2k y Re{44*e 2 ^ +k ^)}e~ 2k ^ 

= ^{T(|4r + KIV^ + AV 2 + 2 |4| |4| k v cos(2£ ± + 2k r x + 0± - 0±)} 

where 0^ = Arg(pf) (p = a,d). Above expressions contain a wealth of information. It 
is not difficult to see that for any values of coefficients 4 an d 4 the fi rs t term always 
dominates so the currents of the two valleys are not only in opposite directions, their direc- 
tions are fixed and do not change with sign of k y . This is unusual and unlike the case of 
ordinary systems described by usual Schrodinger equation where in case of invariance along 
y-direction y-component of current Jy always follow the sign of the k y and can always be 
written as Jy = r ^-p(x) where p{x) is probability density and other factors have obvious 
meanings. Above result shows that, on one hand we can use these states for valley depen- 
dant functionalities by applying a voltage difference along the interface and on other hand, 
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coupled with the fact that current has to be continuous at the interface it implies that the 
states of the two valleys in CR propagating towards the interface at any positive or negative 
angle will carry currents in the same fixed directions close to the interface on its both sides. 
Far from the interface in CR flux have to follow the sign of k y . The strip-like region around 
the interface where valley dependant properties show up are defined by localization lengths 
of the evanescent states on both sides. This also shows that the localized states in IR at 
E inside the band gap modify the behavior of those in CR that usually do not distinguish 
between the two valleys and currents carried by them follow the sign of k y . Another in- 
teresting point is that J^ y is non-zero even when k y = and more important for practical 
purposes is that the factors C ± switch sign with the polarity of the gates. This is clear 
from their expressions where all terms on right side are positive definite except (U — E) 
which is positive for U > and negative otherwise. Thus we can control the directions of 
the currents of the two valleys with the polarity of the gates. This broadens the possible 
set of functionalities we can achieve using these states. Further, since all states with energy 
inside the band gap have similar behavior, smearing of fermi surface due to temperature 
or any intravalley scatterings, which are almost unavoidable in real systems, is harmless 
so experimental realization of a device where valley based operations could be performed 
is expected to be very simple and easy using these localized states. One way is, similar 
to ref 9j as mentioned earlier, to use just the narrow strip region of the interface where a 
small voltage difference along it can be used to get a net valley polarized current or filter 
a valley which can be of particles belonging to either valley depending on the polarity of 
gates used to create the band gap. Two such set ups in series can be used as a valley valve 
where particles can only pass if the polarities of the gates are the same for the two. Another 
way is to use the CR as the source of particles and again in this case depending on the 
polarity of the gates one of the contacts on the lateral sides will collect only particles of one 
valley. These contacts do not need to be very fine on the IR, however, they must not extend 
in CR more than the localized states otherwise they will also collect particles of the other 
valley. And similar to the above case a combinations of such set ups can be used for other 
functionalities. 

For completeness, let's compute the currents J^ y carried by the localized states on IR 
as well as the currents along the interface in CR when the particles in CR travel to- 
wards the interface at a positive energy E and with y-component of wavevector k y . The 
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wavefunctions in CR in this case can be written as ^f(x,y) = ^f(x)e lkyy with ^f(x) = 

e -ikx _|_ c ± | J 6 +kx yjliQYQ ^± an( j c ± are com pl ex coeffi- 

, Vv 

cients, = ( y/l + sin 2 0^sin 0) 2 = fc = ^/E - fcjj , « = y/E + fc* and = tan -1 (^). 

Using the matching conditions (x = 0) = ^(a? = 0) and d x ^f(x)\ x= Q = d x ^f(x)\ x= o 
we determine all the coefficients bf, cf , af and df . As expected for energies inside the band 
gap, \bf | = 1 so there is no net current along x-direction. Further, we find that \af\ = \df\ 
= la^l so we can write jf y = 2 h] ^ ] {U ~§ ± R± |a ± | 2 {=F^/l + (k r /k y ) 2 + sgn(k y ) cos(29± + 2k r x + 
0a — 4>d)}e~ 2kiX . This expression clearly shows the independence of directions of Jf y from 
the sign of k y and other features mentioned above. Currents of the two valleys in CR jf y are 
determined by using tyf (x, y) with above calculated coefficients bf and cf. Figure(l) shows 
the contour plots of the currents of the two valleys along the y-direction J^, where equals 
jf y for x < and J^ y otherwise, for E = 17meV and U = ±50meV as a function of k y 
and distance from the interface. Purple (dark) and off-white (bright ) colors show currents 
along — y and +y respectively. Figures(la,b) show J y for U = ±50meV and U = ^f50meV, 
also indicating that is the same for a given polarity of gates as J y for the opposite 
polarity and this relation holds at all positions x and for all values of k y . It is clear that 
for U = +50meV, close to the interface, J y flows towards — y for all possible positive and 
negative values of k y where as J y flows towards +y for all possible positive and negative 
values of k y . Further, for the opposite polarity of gates, i.e., for U = — 50meV, currents 
of both valleys switch directions. From the figure, the width of the strip-like region where 
this behavior is observed is approximately 15nm for the parameters used so for an interface 
sharp enough that these localized states can be used but smooth enough that intervalley 
scatterings can be ignored is clearly possible. Far from the interface currents for the two 
valleys retain the usual behavior. In IR they vanish when kiX » 1 where as in CR for 
k\x\ » 1 their directions follow the sign of k y . We can also see this fact by calculating jf y 
far from the interface where the localized states have negligible effects. They are given by 
jf = -^-{1 — cos(2A;a; — 20 — <pf)} where 0^ = Arg(bf) which in general may be different 
for the two valleys. These expressions clearly show that far from the interface in CR currents 
of both valleys have the same direction determined by the sign of k y . In figure (2c) sum of 
the two currents is plotted which shows that angular symmetric incidence of the particles 
of the two valleys will result in no net charge current just like the case of common systems, 
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however, of course, there will be two non-zero valley currents flowing in opposite directions. 

Let's summarize above findings and draw some conclusions. We see that the wavefunction 
matching conditions remain the same in case of energy dependant effective mass in gaped 
bilayer graphene system which is very unusual. We showed that the states localized near an 
interface between conducting and insulating regions of BLG on the insulating side despite 
having plane wave parts along the decay direction carry no currents along it. Further, we 
have seen that their properties in transverse direction are strikingly different than those of 
evanescent waves in ordinary systems where direction of current follow the sign of corre- 
sponding wavevector component. Moreover, the states belonging to the two valleys show 
contrasting behavior and their presence also modifies the properties of the localized states 
on the conducting side of the interface. Their valley dependent properties are insensitive 
to various mechanisms that can possibly pose problems for earlier proposals to obtain basic 
valley based functionalities. They are easy to exploit in realistic situations to obtain valley 
polarization, filter particles of a desired valley or for valley switching without dealing with 
the propagation angles of the particles or facing any other experimental challenges. Thus we 
expect them capable of playing an important role in valleytronics based on bilayer graphene 
system. 
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Figure captions: 

Fig.l(color online). Variation of currents Jy of the two valleys K(+) and K'(— ) with 
distance from the interface x and component of wavevector along the interface k y . Figure 
shows that for a given polarity of gates Jy and J~ have fixed opposite directions that do not 
depend on the sign of k y . (a) for U = ±50meV (b) for U = ±50meV and (c) sum 
of the two currents Jy, Jy + Jy for U = ±50meV. From the figure the width of the region 
that supports the valley polarized currents along the fixed directions is roughly 15nm. 
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